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ABSTRACT 

We discuss bounds on the strength of the magnetic fields that could be buried in 
the radiative zone of the Sun. The field profiles and decay times are computed for all 
axisymmetric toroidal Ohmic decay eigenmodes with lifetimes exceeding the age of the 
Sun. The measurements of the solar oblateness yield a bound < 7 MG on the strength 
of the field. A comparable bound is expected to come from the analysis of the splitting 
of the solar oscillation frequencies. The theoretical analysis of the double diffusive 
instability also yields a similar bound. The oblateness measurements at their present 
level of sensitivity are therefore not expected to measure a toroidal field contribution. 

Subject headings: Sun: magnetic fields — Sun: interior 


1. Introduction 

It has been speculated that the solar radiative interior may contain a buried magnetic field, 
which does not penetrate into the convective zone and hence is very difficult to observe. Such a 
field could have formed, for instance, during the early stages of the solar evolution, perhaps as a 
result of the differential rotation which could have stretched the lines of a weak primordial seed 
field (Dicke 1979). Once formed, the large scale field structures would survive until the present 
time because of the high conductivity of the plasma in the radiative zone. 

The observation that the lifetime of the large scale field eigenmodes in the solar interior may 
exceed the solar age was first made by Cowling over half a century ago (Cowling 1945). Cowling’s 
conclusion was confirmed by other authors, in particular by Wrubel (1952) who computed the 
lifetimes of the three longest living poloidal modes and by Bahcall & Ulrich (1971) who computed 
the lifetime of the lowest mode using an early solar model. 
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The motivation for discussing the hypothetical radiative zone (RZ) field changed with time. 
The RZ fields, with a characteristic strength of 10 8 G, were invoked by several authors (Bartenwerfer 
1973; Chitre, Ezer, & Stothers 1973) in the 1970’s as a possible explanation for the deficit of solar 
neutrinos reported by the Homestake experiment. The fields of comparable strength was postulated 
by Dicke (Dicke 1978, 1979, 1982) to explain the so-called Princeton solar oblateness measurements. 
Recently, Gough & McIntyre (1998) argued that a relatively weak, but nonzero, RZ poloidal field 
is necessary to explain the uniform rotation of the radiation zone. Separately, the fields of 30 
MG strength were proposed as a way to improve the agreement between the solar model and the 
measurements of the 8 B neutrino flux at the SNO and SuperKamiokande experiments (Couvidat, 
Turck-Chieze, & Kosovichev 2002). In a still separate development, it was found that the presence 
of a submegagauss toroidal field in the radiative zone could explain all presently available solar 
neutrino data, if neutrino possesses a relatively large transition magnetic moment (Friedland & 
Gruzinov 2002). 

Regardless of the motivation, it is important to establish what bounds, both observational and 
theoretical, exist on the strength of the radiative zone fields. This is the goal of the present paper. 
Both Cowling and Wrubel in their analyses assumed a poloidal topology for the primordial field. 
Such field cannot be stronger than a few gauss at the top of the radiative zone, otherwise it would 
penetrate in the convective zone (CZ) and cause a polarity asymmetry between the two halves of 
the magnetic cycle (Levy & Boyer 1982; Boyer Sz Levy 1984; Boruta 1996). The absence of the 
asymmetry places a very strong bound (B < 10 2 — 10 3 G) on the strength of the poloidal field 
allowed in the solar core. At the same time, the field which is entirely confined to the radiative 
zone, in particular a toroidal axisymmetric field, is not subject to this bound, and may in fact be 
very difficult to observe. 

Before presenting the constraints on the strength of the toroidal field, in Section 2 we compute 
the Ohmic decay eigenmodes for the axisymmetric toroidal field configurations and the correspond¬ 
ing lifetimes, and give the profiles and lifetimes for all modes whose lifetimes exceed the solar age. 
To the best of our knowledge, such a calculation has not been done previously. The problem was 
touched on by Dicke (1982), who presented the profile of the l = 2 mode (which is not the longest 
living mode), without specifying its lifetime. The discussion of the bounds on the field strength 
follows in Section 3. This is followed in Section 4 by the discussion of theoretical instabilities that 
may provide bounds on the strength of the field. In Section 5 we summarize our conclusions. 


2. Eigenmodes 

We will approximate the electric conductivity in the radiative zone of the Sun by the Spitzer 
formula (Spitzer 1962) 

2 5 / 2 (kTf/ 2 
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where m e and e are the electron’s mass and charge, Z e g = is the effective ion 

charge, T is the temperature, 7 E is an order one coefficient describing the deviation of the plasma 
from the ideal Lorentz gas, and A is the Coulomb logarithm. Numerically, this yields 

<7 = (2.37 x 10 17 s -1 )Tg /2 —(2) 

Z e ffX 

where Tq is the temperature in millions of Kelvins. The quantity Z e Q can be computed by using the 
tables of the solar model; its dependence on the solar radius for the BP2000 solar model (Bahcall, 
Pinsonneault, & Basu 2001) is shown in Fig. 1(a). The Coulomb logarithm A for Tq > 0.6 is given 
by 

A = 2.4 — 1.15 log 10 n + 2.3 log 10 T e , (3) 

where the number density of electrons n is in units of Avogadro number per cm 3 . Finally, the 
coefficient 7 E for Z e g- in the interval 1 < Z e g < 2 is well described by 

1e = 0.58 + 0.33 log 10 Z eff . (4) 

The radial dependence of the factor ^e/Z^X in the Sun is shown in Fig. 1(b). 

The principal approximation affecting the accuracy of Eq. (2) is using the Spitzer conductivity 
in the center of the Sun, where Coulomb logarithm is not large, A ~ 2.9. 

The corresponding magnetic diffusivity, 7 = c 2 /(47 tct) is given by 

rj = (3.0 x 10 2 cm 2 s -1 )^5—T 6 -3/2 . (5) 

7 E 

In the convective zone, the effective magnetic diffusivity is very large because of turbulence. 

The Ohmic decay of the field is governed by the magnetic diffusion equation, 

<9(B = —V x (r/V x B), (6) 

which follows from <9fB/c = —V x E, E = j/cr, and V x B = 47rj/c. Eqs. (5,6) can be used to 
estimate the minimal size Z m ; n of the magnetic field features that survive over the age of the Sun, 
Z m in ~ vW, where t & = 4.6 x 10 9 years is the age of the Sun. Since ?? varies with the distance from 
the center, so does / min - At r = 0.6 Rq one finds Z m j n ~ R@/ 9, while at r = 0.25 Rq Z m ; n ~ Rq/ 20. 
Thus, field modes with the size comparable to the size of the radiative zone are expected to survive 
until the present time. 

To make the preceding rough estimate more precise, we next compute the profiles and the 
lifetimes of the axisymmetric toroidal eigenmodes that are localized to the solar RZ. For an ax- 
isymmetric toroidal field B = B(r, 0,t)e^, Eq. (6) gives 


d t B = r](v 2 B - 1 B) + r/r 1 d r (rB ) 

\ r z sm 0 J 


( 7 ) 
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where ' = d/dr. The eigenmodes have the form B = e t ^ Tn > l F n i(r)P/(cos 9) where n = 1,2,..., 
I = 1,2,.... The function F n j satisfies the equation 

r] (F" + 2 r~ l F' - 1(1 + 1 )r~ 2 F) + rj'(F' + r~ l F ) = -r z _1 F. (8) 

The boundary conditions are F = 0 at r = 0 and r = i?cz- 

The lifetimes of various modes, computed for the BP2000 Sun (Bahcall, Pinsonneault, & Basu 
2001), are tabulated in Table 1. The n = 1, / = 1 mode has the longest lifetime, rip = 24Gyr. As 
can be seen from the table, there are eight modes whose lifetimes exceed the age of the Sun and 
three modes whose lifetimes exceed 10 billion years. Therefore, the toroidal field in the radiation 
zone of the Sun can in principle have complex structure. The profiles F n j(r) of the eight longest 
living modes are shown in Figs. 2 and 3 and tabulated in Table 2. 

Although the preceding calculation does not take into account the changes of the solar pa¬ 
rameters with time, the effect of the solar evolution on the lifetimes is expected to be small. As a 
result of hydrogen burning, as the Sun evolves from t = 0 to t = 4.7 Gyr, the core contracts and 
becomes hotter, while the outer layers of the RZ expand and cool. Since the lifetime depends on 
the product Z 2 T 3 / 2 , there is partial cancellation, both in the core and in the outer part of the RZ, 
between the effects of changing temperature and changing size. Moreover, the crossover radius, at 
which the temperature and the radius distance do not change, is at about O.15i?0 (Demarque & 
Guenther 1991), close to where the maximum of the n = 1, l = 1 mode. 


3. Observational upper bounds on the magnetic field 

In this Section we discuss observational upper bounds on the allowed strength of the magnetic 
field in the radiation zone. A variety of independent arguments rule out fields in excess of < 10 s 
G. We briefly summarize the arguments from helioseismological measurements of the sound speed 
and from the measurements of the 8 B neutrino flux. A much stronger bound on the axisymmetrical 
fields comes from the measured oblateness of the Sun. We derive a bound on the strength of the 
held specializing to the longest living n = 1, l = 1 mode. We also point out that a comparable 
bound could be obtained from the helioseismological measurements of the mode splitting. 


3.1. Bounds from helioseismology and S B neutrino flux 

The magnetic held of strength B adds a contribution to pressure p m = B 2 /8ir, correspondingly 
decreasing the gas pressure p g . By estimating the change in the gas temperature as | ST/T\ rsj Pm/Pg, 
one obtains a bound on the allowed strength of the held in the core and in the RZ from the observed 
hux of the 8 B neutrinos and the helioseismological measurements of the sound speed. The two 
constraints are complimentary, because the 8 B neutrinos are produced very close to the center of 
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the Sun (r < 0.09 Rq), where the accuracy of the helioseismological measurements is relatively 
poor. 


The flux of the 8 B neutrinos neutrinos is not directly tied to the solar luminosity and has a 
strong dependence on the core temperature, c/>( 8 B) oc T@, where /3 ~ 24 (Bahcall & Ulmer 1996). 
Therefore, the flux </>( 8 B) can be used as a probe of the central temperature. The flux </>( s B) has 
been measured by the SNO experiment using a neutral current reaction is + d —* is + p + n and 
assuming oscillations into active neutrinos. The measured value, 5.09 ± 0.64 x 10 6 cm 2 s 1 , is in 
good agreement with the predictions of the standard solar model, 5.05+g'g x 10 6 cm -2 s^ 1 (Bahcall, 
Pinsonneault, & Basu 2001). Combining the above ingredients, we obtain an order of magnitude 
bound on the strength of the fields in the core, 



( \ 1/2 
V <K 8 B) ) 


~2x 10 8 G. 


( 9 ) 


A similar order of magnitude bound follows from the measurements of the sound speed. Taking 
the uncertainty on the sound speed to be <5c s /c s 10 3 we find for r = 0.2 Rq, 




~ 0.4 x 10 8 G. 


( 10 ) 


To obtain more accurate bounds in both cases requires incorporating magnetic fields in the 
solar model calculations. Recent analysis by Couvidat, Turck-Chieze, & Kosovichev (2002) finds 
the limit B < 3 x 10' G from the measurements of the sound speed. 


3.2. Constraints from Solar Oblateness 

It is well known that the rotation distorts the shape of the Sun, making it oblate. The 
theoretical value of the solar rotational oblateness e = 1 — R po ie/-^equator is e rot = 9 x 10“ 6 (Godier 
& Rozelot 2000). The magnetic fields will further distort the equilibrium figure of the Sun. If we 
take the recently measured oblateness of the Sun to be (10 ± 3) x 10 -6 (Godier & Rozelot 2000), 
we must require that the contributions to the oblateness from the magnetic fields not exceed the 
purely rotational oblateness by over ~ 30%, which translates into a bound on R max . 

The distortion is obtained from the equilibrium equation 

V(p + p m ) + %£ + pV^ = 0. (11) 

Here p is the gas pressure, p rn = B 2 /(8n) is the pressure of the toroidal magnetic field, R = r sin# 
is the cylindrical radius, R is the unit vector along cylindrical radius, p is the density, (j) is tl ie 
gravitational potential. We can exclude p by taking the curl, and then linearize, that is leave only 
first order terms in p m : 


4>'q d e Sp - p' 0 d e 5(j) = 2 cot 9d r p m - 2 r 1 d e p m . 


( 12 ) 
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Here 5p and S(p are density and potential perturbations due to magnetic pressure, po and (j )o are 
the unperturbed density and potential of the Sun, and prime is the derivative with respect to r. 

We now specialize to the longest living mode, p m (r, 9) = B 2 nax [F\] (r)] 2 sin 2 9/8n , with Fi i(r) 
tabulated in Table 2. Then dp and 5(j> are oc T^cos 9) and, using the Poisson equation V 2 (j) = AnGp, 
we obtain 

(47rG)-V{,(<ty" + 2r~ V - 6r" 2 <ty) - p' o 50 = -(4/3)bL*^i,i(*i,i - r^F^). (13) 

This equation was solved numerically, and the calculated oblateness of the Sun was found to be 

emagn = -6.4 X 10- 8 . (14) 

By requiring that e magn < 3 x 10 —6 , we get an upper bound 

£ max < 7MG. (15) 

Thus the bound from oblateness is much stronger then the ones considered earlier. 

The effect of the 12 MG field would be to completely cancel out the solar rotational oblateness, 
making the Sun into a perfect sphere. The fields of the strength 3 x 10' G, that were invoked by 
Couvidat, Turck-Chieze, & Kosovichev (2002) as a possible way to improve the agreement between 
the observed and predicted S B fluxes, would make the figure of the Sun strongly prolate and hence 
are not allowed. 


3.3. Splitting of Solar Oscillation Frequencies 


The observations of the splitting of solar oscillation frequencies can provide another way to 
constrain the strength of the magnetic field. Antia, Chitre, & Thompson (2000) analyzed the 
observed splittings and derived an upper bound of B < 0.3MG for a toroidal field of radial extent 
d = O.O4f?0 near the base of the convective zone at r = 0.7i?g. We can write their upper bound in 
a form 1 


d B 2 
r 8irp(r) 


< 3 x 10“ 6 . 


(16) 


If we assume that the error bars on the measured splitting coefficients do not increase too much as 
we move deeper into the Sun, we obtain an upper bound on the strength of the lowest field mode 
that is comparable to the one derived from oblateness. To make this bound more precise requires 
repeating the full analysis of Antia, Chitre, & Thompson (2000) with the magnetic field profile 
composed of the modes found in Sect. 2. 


lr This is a much stronger upper bound than the upper bound B 2 /(8ivp) < 10 3 , considered in Sect. 3.1. 
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4. Theoretical Considerations 

In this section we describe theoretical bounds on the magnitude of the radiative zone magnetic 
field. The bounds are based on the observation that a sufficiently strong field will rise to the surface 
as a result of thermal diffusion in the solar plasma. We first consider the simplified scenario in 
which the field rises as a whole and then address the possibility that the field may fragment and the 
fragmented flux tubes rise to the surface. The latter analysis gives a particularly stringent bound 
on the field strength. We also comment on the role of the ideal magnetohydro dynamic instability 
that may exist for the field configuration in question. 


4.1. Constraint from Magnetic Buoyancy 


A robust bound on the strength of the primordial magnetic field follows from the consideration 
of the rise time in the stably stratified RZ (Parker 1974, 1979). The essence of the effect is that 
a flux tube in hydrostatic equilibrium in a stably stratified medium has a lower temperature than 
the temperature of the surrounding medium. The resulting flow of heat into the tube causes it to 
rise to assume a new equilibrium position. The rise velocity u thus depends on the field strength 
(through temperature gradient) and the heat conduction properties of the plasma, 


IB 2 / A \ 2 / 
6 8irp \RJ p 


(17) 


Here p is the equilibrium pressure, I = Lq/A'kt 2 is the heat flux, R is the radius of the flux tube, 
A = d In T/dr is the temperature scale height of the medium, and 5 = ( dlnn/dlnT ) — 1/(7 — 1) is 
a parameter determining the stability of the medium (for an adiabatic medium 5 = 0). 

Numerically, the rise time from r = 0.2 R & to the top of the RZ is 


f rise ~ 3 x 10 7 yr 



(18) 


so that the field B ~ 10 s G would rise in the time comparable to the age of the Sun. For the field 
B ~ 0.5 x 10 6 G considered in Friedland & Gruzinov (2002) the distance traveled over the age of 
the Sun is only about 8 km and the effect can be safely ignored. 

The above mechanism can be used to make an argument that the primordial RZ field could 
not have formed with the strength of 10 8 G. This is possible using the abundance of Beryllium at 
the solar surface. Be is destroyed at T = 3 x 10 6 K, which is the temperature at r = 0.61 -Rq. Since 
Beryllium observed at the solar surface is not noticeably depleted, the material at r ~ 0.7 Rq and 
the material at r ~ 0.6 -Rq have not been mixed during the lifetime of the Sun (except for a brief 
period during the pre-main sequence convection stage). The rise of the 10 8 G field would have 
forced such mixing, and therefore could not have happened. 



4.2. Double Diffusive Instability 


The preceding analysis showed that the rise velocity of a flux tube is inversely proportional to 
the diameter of the tube squared. The simple reason for this is that small flux tubes are heated 
faster than large ones. It is therefore conceivable that instead of rising as a whole, the magnetic 
field might fragment into small flux tubes, which would rise to the surface faster. This instability 
is counteracted by the diffusion of the magnetic field lines. Superficially, it may appear that, since 
the heat diffusivity of the solar plasma, 


16 aT 3 
2>np 2 c v ’ 


(19) 


is some four orders of magnitude greater than the magnetic diffusivity rj, given in Eq. (5), any 
field configuration will be unstable. However, a more careful consideration of the effect shows that 
the relevant parameter controlling the instability is the combination ( x/v)(Pm/Pg )• The stability 
analysis was carried out by Schubert (1968), who found that for the field configuration to be stable, 
it is necessary that 

dlnjpp-'i) > x B 2 9 In (Bj/pr) 
dr rj 47rp g dr 

Substituting in this condition the profile of the longest living toroidal eigenmode (1 = 1, n = 1), we 
find that if the normalization of the field, B max , exceeds ~ 2.1 MG, the field configuration becomes 
unstable at r ~ 0.6 -Rq. Therefore, this argument provides the most stringent bound on the toroidal 
magnetic field. 


4.3. Ideal magnetohydrodynamic instability 

An arbitrarily weak purely toroidal magnetic field is unstable to ideal magnetohydrodynamic 
(MHD) perturbations (Tayler 1973; Goossens & Veugelen 1978; Goossens & Tayler 1980). How¬ 
ever, one can show (Gruzinov & Friedland 2002) that the perturbation considered by Goossens & 
Veugelen (1978) is stabilized by a poloidal field B p which is much weaker than the destabilizing 
toroidal field B t {B 2 /8 ttp < (B 2 /8irpj 2 ). 

The general problem of ideal MHD stability of toroidal plus poloidal magnetic field in spherical 
geometry with stable stratification is unsolved. One also notes that toroidal fields satisfying obser¬ 
vational upper bounds are very weak ( B 2 /(8irp ) <3x ICE 6 ), and the unstable perturbations have 
a very small wavelength (Goossens & Veugelen 1978). Then even a trace of differential rotation in 
the RZ can affect stability in an unknown manner. It is therefore difficult to extract a concrete 
bound on the strength of the magnetic field from the analysis of the MHD instability with the 
present level of understanding of the problem. 
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5. Discussion 

We have discussed a variety of constraints on the allowed strength of the magnetic fields in the 
solar radiative zone. We have shown that the strongest constraints come from the measurements 
of the solar oblateness, from the splitting of the oscillation frequencies, and from the theoretical 
analysis of the double diffusive instability of the field. All three constraints point to maximal field 
values of less than a few magegauss. That very similar bounds emerge from the analysis of such 
physically different phenomena is in itself quite remarkable. Alternatively, our results show that 
magnetic field that would measurably distort the figure of the Sun would cause an instability in 
the Sun, and hence is not allowed. This explains why the oblateness measurements at their present 
level of sensitivity have not seen a toroidal field contribution. 

Our analysis excludes the possibility of perturbing the solar model with the fields of the order 
of several tens of megagauss that have been invoked to obtain a better agreement with the solar 
neutrino data (Couvidat, Turck-Chieze, & Kosovichev 2002). 

The bounds obtained here have a direct implication for the analysis of the neutrino survival 
probability discussed in Friedland & Gruzinov (2002). The survival probability depends on the 
value of the product of the neutrino transition moment, /j, and the magnetic field normalization, 
Umax' The value of the product, which was found to give a good fit to the solar neutrino data, is 
/i-£> max f ~ 0.5 x HT 11 MGxgB (/i b is the Bohr magneton). The upper bound on the magnetic field 
strength therefore yields a lower bound on the size of the transition moment that is necessary to 
fit the solar neutrino data. The bound from oblateness yields /i > 7 x lCK 13 ^# , while the limit 
from the double diffusive instability gives /i > 2 x lCP 12 ^#. 

Among the theoretical issues that require further study is the problem of ideal MHD stability 
described in Section 5 and the possibility that the symmetry axis of the magnetic field may become 
tilted as a result of the dissipative motions of the plasma (Spitzer 1958; Mestel, Nittmann, Wood, 
& Wright 1981). 

In summary, the bounds on the strength of the hypothetical magnetic fields in the radiative 
zone of the Sun described here may be useful for searches for such fields and also for the discussions 
of their physical implications. 


We thank John Bahcall for valuable discussions and his encouragements during the course of 
this work. We are also grateful to Axel Brandenburg, Sarbani Basu, and Peter Goldreich for useful 
discussions. A. F. was supported by the W. M. Keck Foundation. 
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Fig. 1. The radial dependence of the effective ion charge Z e Q (left) and the quantity ^e/Z c q\ 
(right), computed for the BP2000 Sun. 
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Fig. 2.— The radial dependence F Uj i(r) of the l = 1 eigenmodes (left) and l = 2 eigenmodes (right) 
whose lifetimes exceed the age of the Sun. 



Fig. 3.— The radial dependence F n j(r) of the l = 3 (left) and l = 4 (right) eigenmodes whose 
lifetimes exceed the age of the Sun. 
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1 

2 

3 

4 

5 

1 

24.3 

12.8 

7.7 

5.2 

3.7 

2 

10.0 

6.6 

4.6 

3.4 

2.6 

3 

5.3 

4.0 

3.0 

2.4 

1.9 

4 

3.3 

2.6 

2.1 

1.7 

1.4 

5 

2.2 

1.8 

1.5 

1.3 

1.1 


Table 1: The lifetimes of the field eigenmodes in billions of years as a function of n (vertical) and l 
(horizontal). 
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Table 2. 


The radial dependence T) n (r) of the eigenmodes whose lifetimes exceed the age of the 

Sun. 


r 

Fi, i 

F2,l 

7*3, i 

Fi,2 

7*2,2 

7*1,3 

7*2,3 

7*1,4 

0.01 

0.1130 

0.1669 

0.2239 

0.0077 

0.0147 

0.0005 

0.0011 

0.0000 

0.03 

0.3326 

0.4809 

0.6250 

0.0672 

0.1271 

0.0120 

0.0277 

0.0020 

0.05 

0.5325 

0.7378 

0.8978 

0.1772 

0.3228 

0.0523 

0.1164 

0.0146 

0.07 

0.7019 

0.9118 

0.9994 

0.3217 

0.5537 

0.1306 

0.2759 

0.0501 

0.09 

0.8341 

0.9929 

0.9345 

0.4810 

0.7673 

0.2455 

0.4839 

0.1184 

0.11 

0.9265 

0.9862 

0.7453 

0.6363 

0.9219 

0.3863 

0.6978 

0.2216 

0.13 

0.9804 

0.9081 

0.4923 

0.7723 

0.9948 

0.5372 

0.8728 

0.3528 

0.15 

0.9997 

0.7809 

0.2339 

0.8790 

0.9834 

0.6816 

0.9769 

0.4987 

0.17 

0.9902 

0.6270 

0.0129 

0.9519 

0.9006 

0.8061 

0.9972 

0.6435 

0.19 

0.9581 

0.4661 

-0.1477 

0.9910 

0.7671 

0.9017 

0.9390 

0.7727 

0.21 

0.9096 

0.3128 

-0.2429 

0.9995 

0.6062 

0.9646 

0.8200 

0.8761 

0.23 

0.8504 

0.1766 

-0.2799 

0.9824 

0.4381 

0.9955 

0.6631 

0.9486 

0.25 

0.7850 

0.0623 

-0.2722 

0.9456 

0.2782 

0.9977 

0.4907 

0.9890 

0.27 

0.7171 

-0.0288 

-0.2348 

0.8944 

0.1365 

0.9760 

0.3206 

0.9999 

0.29 

0.6493 

-0.0976 

-0.1810 

0.8339 

0.0179 

0.9358 

0.1655 

0.9852 

0.31 

0.5834 

-0.1465 

-0.1216 

0.7679 

-0.0759 

0.8823 

0.0329 

0.9500 

0.33 

0.5207 

-0.1785 

-0.0640 

0.6997 

-0.1460 

0.8198 

-0.0742 

0.8995 

0.35 

0.4618 

-0.1967 

-0.0128 

0.6316 

-0.1946 

0.7521 

-0.1556 

0.8383 

0.37 

0.4072 

-0.2038 

0.0295 

0.5652 

-0.2250 

0.6823 

-0.2130 

0.7703 

0.39 

0.3570 

-0.2026 

0.0621 

0.5016 

-0.2402 

0.6125 

-0.2493 

0.6990 

0.41 

0.3111 

-0.1951 

0.0852 

0.4417 

-0.2434 

0.5444 

-0.2679 

0.6268 

0.43 

0.2695 

-0.1832 

0.0996 

0.3859 

-0.2375 

0.4794 

-0.2723 

0.5559 

0.45 

0.2318 

-0.1685 

0.1067 

0.3344 

-0.2248 

0.4180 

-0.2655 

0.4875 

0.47 

0.1979 

-0.1520 

0.1076 

0.2873 

-0.2074 

0.3610 

-0.2507 

0.4229 

0.49 

0.1676 

-0.1347 

0.1039 

0.2445 

-0.1871 

0.3084 

-0.2301 

0.3627 

0.51 

0.1406 

-0.1173 

0.0967 

0.2058 

-0.1653 

0.2605 

-0.2061 

0.3073 

0.53 

0.1165 

-0.1003 

0.0872 

0.1712 

-0.1430 

0.2172 

-0.1801 

0.2567 

0.55 

0.0953 

-0.0840 

0.0762 

0.1403 

-0.1209 

0.1783 

-0.1536 

0.2111 

0.57 

0.0765 

-0.0689 

0.0645 

0.1128 

-0.0998 

0.1436 

-0.1275 

0.1702 

0.59 

0.0599 

-0.0549 

0.0527 

0.0886 

-0.0799 

0.1128 

-0.1026 

0.1339 

0.61 

0.0454 

-0.0421 

0.0412 

0.0672 

-0.0616 

0.0856 

-0.0794 

0.1017 

0.63 

0.0327 

-0.0306 

0.0304 

0.0484 

-0.0449 

0.0617 

-0.0580 

0.0734 

0.65 

0.0216 

-0.0204 

0.0204 

0.0320 

-0.0299 

0.0408 

-0.0388 

0.0487 

0.67 

0.0120 

-0.0113 

0.0114 

0.0178 

-0.0167 

0.0227 

-0.0217 

0.0272 

0.69 

0.0037 

-0.0035 

0.0035 

0.0055 

-0.0051 

0.0070 

-0.0067 

0.0087 
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